1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


Topographic and stochastic influences on pahoehoe lava lobe emplacement 


Christopher W. Hamilton 1 ' 2 , Lori S. Glaze 1 , Mike R. James 3 , and Stephen M. Baloga 4 
1 Planetary Geodynamics Laboratory, NASA Goddard Space Flight Center, Greenbelt, 
MD, 20771, USA, Christopher .hamilton@nasa. gov . 

“Department of Astronomy, University of Maryland, College Park, MD, USA. 

Lancaster Environment Centre, Lancaster University, Lancaster, UK. 

4 Proxemy Research, Laytonsville, MD, USA 

Abstract 

A detailed understanding of pahoehoe emplacement is necessary for developing 
accurate models of flow field development, assessing hazards, and interpreting the 
significance of lava morphology on Earth and other planetary surfaces. Active pahoehoe 
lobes on Kllauea Volcano, Hawai‘i, were examined on 21-26 February 2006 using 
oblique time-series stereo-photogrammetry and differential global positioning system 
(DGPS) measurements. During this time, the local discharge rate for peripheral lava lobes 
was generally constant at 0.0061 ± 0.0019 m /s, but the areal coverage rate of the lobes 
exhibited a periodic increase every 4.13 ± 0.64 minutes. This periodicity is attributed to 
the time required for the pressure within the liquid lava core to exceed the cooling- 
induced strength of its margins. The pahoehoe flow advanced through a series of down- 
slope and cross-slope breakouts, which began as ~0.2 m-thick units (i.e., toes) that 
coalesced and inflated to become approximately meter-thick lobes. The lobes were 
thickest above the lowest points of the initial topography and above shallow to reverse- 


24 facing slopes, defined relative to the local flow direction. The flow path was typically 

25 controlled by high-standing topography, with the zone directly adjacent to the final lobe 

26 margin having an average relief that was a few centimeters higher than the lava- 

27 inundated region. This suggests that toe-scale topography can, at least temporarily, exert 

28 strong controls on pahoehoe flow paths by impeding the lateral spreading of the lobe. 

29 Observed cycles of enhanced areal spreading and inflated lobe morphology are also 

30 explored using a model that considers the statistical likelihood of sequential breakouts 

3 1 from active flow margins and the effects of topographic barriers. 

32 

33 1 . Introduction 

34 Basalt is the most common rock type on the surface of terrestrial bodies 

35 throughout the solar system and — by total volume and areal coverage — pahoehoe flows 

36 are the most abundant form of basaltic lava in subaerial and submarine environments on 

37 Earth (Self et al., 1994, 1998). Pahoehoe flow fields are composed of flows, lobes, and 

38 toes, which represent a continuum of lava emplacement scales, with toes being the 

39 smallest elements. Pahoehoe flows commonly advance as toes break out along lobe 

40 margins (Hon et al., 1994; Crown and Baloga, 1999; Hoblitt et al., 2012). These new toes 

4 1 quickly cool and develop a rheological gradient that includes an inferred three-part 

42 structure composed of a brittle outer crust, underlying viscoelastic layer, and inner molten 

43 core (Hon et al., 1994). As toes develop sufficient strength to retain incoming lava, they 

44 can pressurize, coalesce, and inflate to form lobes that are interconnected with other 

45 portions of the flow through internal fluid pathways (Walker, 1991, 2009). The evolution 

46 of these pathways can vary considerably (Kauahikaua et al., 1998), but the initial 


2 



47 emplacement of pahoehoe flow fronts can influence the subsequent development of the 

48 flow. Consequently, a detailed understanding of processes operating on the scale of toes 

49 and lobes can provide important information for modeling aspects of flow field 

50 development, assessing hazards, and interpreting the significance of lava flow 

5 1 morphology on Earth and other planetary bodies (Peitersen and Crown, 2000; Byrnes and 

52 Crown, 2001). 

53 Hon et al. (1994) suggest that “microtopography” (i.e., tens of centimeters relief), 

54 herein termed “toe-scale topography”, can strongly affect the emplacement of pahoehoe 

55 lava flows that reach a final inflated thickness of several meters. Previous studies have 

56 attempted to determine the relationship between topography and pahoehoe emplacement 

57 (Crown and Balgoa, 1999; Peitersen and Crown, 2000; Byrnes and Crown, 2001), but 

58 these studies have been unable to quantify the effects of toe-scale topography at the flow 

59 front due to data resolution limitations. Here, oblique digital photogrammetry and 

60 Differential Global Positioning System (DGPS) measurements are combined to quantify 

61 key physical parameters involved in pahoehoe flow advance through development of toes 

62 and lobes, which constitute the fundamental building blocks of all pahoehoe flows 

63 (Walker, 1991; Self et al., 1996, 1998; Thordarson and Self, 1998). These field 

64 observations are also used to inform and constrain a new stochastic model that describes 

65 pahoehoe lobe emplacement in terms of a balance between random and non-random 

66 processes. 

67 

68 2. Background 
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69 Pahoehoe flow fields exhibit a multitude of branching events at several scales 

70 (Crown and Baloga, 1999) and include self-similar morphologies (Bruno et ah, 1994) that 

7 1 enable local observations to be used to understand larger aspects of flow fields and their 

72 emplacement. However, these patterns are complicated and irregular because of the 

73 influence of stochastic processes operating on multiple scales (Kilburn, 1996). This 

74 makes it important to identify the scales over which recurring processes and structures 

75 combine to form the overall pattern of a pahoehoe flow field. Understanding how 

76 eruption parameters, deterministic factors, and random influences affect the emplacement 

77 and modification of pahoehoe flows, also makes it possible to develop improved 

78 probabilistic models to describe likely patterns of flow growth (Glaze and Baloga, 2013). 

79 This study examines active pahoehoe lava lobes on Kllauea Volcano, Hawai‘i, to 

80 characterize the effects of toe-scale topography on flow emplacement and identify key 

8 1 parameters needed to model statistical aspects of the process. 

82 Pahoehoe flows grow through a combination of areal spreading and inflation, with 

83 new toes breaking out along lobe margins and quickly cooling by radiation (Hon et ah, 

84 1994; Keszthelyi and Denlinger, 1996; Harris et ah, 2007a, 2013). Once a thin skin 

85 develops on the surface of the new toes, cooling becomes increasingly dominated by 

86 conduction (Castruccio et al., 2013). The crust also helps to retain incoming lava, which 

87 slows the advance of the new toes as their internal pressure decreases relative to the 

88 increasing confining strength of the growing crust (Hon et al., 1994). Given a constant 

89 influx of lava, this leads to a reduction in the areal coverage rate of the toes and an 

90 increase in the rate of inflation. Inflation will then continue to dominate the emplacement 

91 process until the internal pressurization of the lobe reaches a threshold required to rupture 
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92 the flow margins and initiate a new cycle of enhanced areal spreading (Hoblitt et ah, 

93 2012). 

94 Within an active pahoehoe flow, internal connections between numerous lobes 

95 can form extensive lava pathways (Kauahikaua et ah, 1998). Preferred pathways within 

96 these lava networks tend to concentrate in existing topographic lows and may develop 

97 into a tube-system as the surrounding flow cools and stagnates. This allows pahoehoe 

98 flows to advance great distances by transporting lava to the flow front through thermally 

99 insulated pathways that can reduce cooling rates to only 0.5-1 .0°C/km (Helz at ah, 

100 1991). As lava is distributed from a master pathway into peripheral lobes and toes, the 

101 local discharge rates will progressively decrease with each branch in the system. Local 

102 discharge rates may therefore be constant over short periods of time (Hon et ah, 1994), 

103 but are expected to evolve as the lobes at the former flow front are incorporated into the 

104 growing body of the lava transport system. In this study, lobes along the periphery of a 

105 low discharge pahoehoe flow are examined to explore flow front emplacement processes 

106 and inflation. 

107 

108 3. Study Area 

109 The Prince Kuhio Kalaniana‘ole (PKK) flow was erupted from Pu‘u ‘0‘o on 

1 10 Kllauea Volcano, Hawai‘i, from March 2004 to June 2007 (Fig. 1). Active margins along 

111 the flow were examined in the field between 2 1-26 February 2006. During this time, lava 

112 from Pu‘u ‘0‘o was primarily transferred though an established tube system that reached 

113 the ocean (Koeppen et al., 2013). However, during February 2006, the master tube- 

1 14 system also fed a series of small surface flows located ~7 km south of Pu‘u ‘0‘o. Harris 
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115 et al. (2007a) referred to this portion of the PKK flow as the “Hook flow”. Field 

116 observations and thermal imagery were used to establish that this flow exhibited 

117 breakouts of S-type pahoehoe (Wilmoth and Walker, 1993) from three small flows on 21 

118 February. However, by 23 February, the activity had dwindled to the distal margin of just 

119 one flow. Waning activity was also evidenced by a decrease in maximum core 

120 temperature from 1 141°C on 21 February, to 1 131°C on 23 February, and a decrease in 

121 areal coverage rate from -0.22 nr/s on February 22 to -0.03 m 2 /s on 23 February (Harris 

122 et al., 2007a). Additionally, there was a corresponding decrease in discharge rate feeding 

123 small surface flows, from -0. 14 m /s on February 22 to -0.03 m /s on 23 February 

124 (Harris et al., 2007a). These surface flows stagnated on February 24. Measurements 

125 reported in this study were made contemporaneously with the observations of Harris et al. 

126 (2007a), who used thermal-image chronometry to establish broad relations among 

127 cooling, discharge rate, and areal coverage rates for the whole flow. 

128 To place the local discharge rates discussed within this study into a broader 

129 context, the Pu‘u ‘O‘o-Kupaianaha eruption began erupting in Kllauea Volcano’s East 

130 Rift Zone (ERZ) in January 3, 1983, and between 1983 and 2002 its time-averaged 

131 discharge rate (Harris et al., 2007b) was -3.81 m /s dense-rock-equivalent lava (Heliker 

132 and Mattox, 2003). However, the Pu‘u ‘0‘o-Kupaianaha eruption has undergone 

133 numerous changes throughout its history (Heliker and Mattox, 2003) and in 2003-2007 

134 the magma supply rate to the ERZ increased significantly during a mantle-derived surge 

135 to Kllauea Volcano (Poland et al., 2012). During this time, the magma supply to Kllauea 

136 Volcano increased by a factor of two (Poland et al., 2012). By 2005 the supply rate to the 

137 ERZ increased to -5.1 1 m /s (Poland et al., 2012) and so at the time of the field campaign 
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138 in February 2006, effusion rates from Pu‘u ‘0‘o were probably higher than the long-term 

139 average. Tilt meter records from the Hawaiian Volcano Observatory’s Pu‘u ‘0‘o cone 

140 station also show that there was an oscillatory pattern of inflation and deflation between 

141 2 1-26 February 2006 (Orr, 2011). This pattern began with an inflation event on 2 1 

142 February, which was followed by gradual deflation over the next five days through a 

143 series of smaller tilt cycles (Tim Orr, personal communication, 2013). Orr (2011) showed 

144 that these tilt cycles were associated with fluctuations in lava discharge through the tube 

145 system and may have contributed to the observed breakouts from the master tube. 

146 Discharge rates associated with the small surface flows examined within this study were 

147 therefore much lower than for the concurrently active PKK tube and the eruption as a 

148 whole. 

149 

150 4. Methods of Digital T errain Modeling 

151 To quantify aspects of the pahoehoe emplacement processes, two complementary 

152 methods were employed. First, time-series oblique stereo-imaging and photogrammetry 

153 (Robson and James, 2007) were used to determine the changing geometric properties of 

154 an active lava lobe, herein referred to as Lobe A. Second, Differential Global Positioning 

155 System (DGPS) measurements of topography before and after the emplacement of a 

156 second pahoehoe lobe (Lobe B) were used to assess the relationship between initial 

157 topography and final lobe morphology. Lobes A and B were located along the same lava 

158 pathway system, with Lobe B located ~10 m directly down-flow of Lobe A. 

159 Time series stereo-image pairs of Lobe A were acquired on 23 February using 

160 two synchronized tripod-mounted digital Single Lens Reflex cameras (6 mega-pixel 
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Canon EOS 300D), which collected image pairs every minute over a duration of 35 
minutes (Robson and James, 2007). The cameras were mounted approximately one meter 
apart and had pre-calibrated, fixed-focus 28 mm lenses. Photogrammetric control was 
provided using spherical targets (0.25 mm in diameter) with positions detennined by 
DGPS measurements. The analysis was carried out using Vision Measurement System 
(Robson and Shortis, http://www.geomsoft.com) with image matching by the GOTCHA 
dense matching algorithm (Gruen, 1985; Day and Muller, 1989; Otto and Chau, 1989). A 
typical stereopair provided topographic point cloud data with positional precisions of 25, 
65 and 20 mm in easting (x), northing (y), and elevation (z), respectively. To quantify 
surface changes, these point cloud data were gridded in Surfer using a Kriging algorithm 
over a 0. 1 m/pixel x-y grid. The data were then imported into ArcGIS. 

For Lobe B, x, y, and z coordinates for topographic surfaces were acquired before 
(23 February) and after lava emplacement. Post-emplacement surveys were performed on 
24-26 February; however, no changes were observed after 24 February and so data from 
these three surveys were combined together. Data were collected at a sampling rate of 0.5 
seconds using two Ashtech Z-Xtreme DGPS receivers — one operating in kinematic mode 
and the other as a static base-station. These data were post-processed using Ashtech 
Solutions 2.60 with differential, L1/L2 (1575.42 MHz/1227.60 MHz) frequency 
correction for ionospheric delay error, and International Global Navigation Satellite 
System Service Standard Product 3 precise orbit file corrections. Resulting data points 
were filtered to remove measurements with more than 0.05 m residual accuracy error. On 
volcanic terrains, Hamilton et al. (2010) estimate the vertical precision of single traverse 
kinematic DGPS surveys to be ±0.08 m, and double traverses being ±0.06 m, with 
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184 horizontal precision being better than the vertical measurements. Most gridded cells in 

185 this dataset contain multiple DGPS observations and should have a vertical precision of 

186 ±0.06 m. The resulting datasets for the initial and final topography include 14,618 and 

187 20,585 records, respectively. To avoid spatial aliasing and to eliminate redundant data, 

188 the BLOCKMODE function in Generic Mapping Tools (GMT; Wessel and Smith, 1991) 

1 89 was used to average z within each pixel. Filtered x, y, and z, data were then interpolated 

190 using NEARNEIGHBOR in GMT to form a 0. 1 m/pixel x-y grid. 

191 

192 5. Time Series Observations of Pahoehoe Lobe Emplacement 

193 Stereogrammetric point clouds and interpolated Digital Terrain Models (DTMs) 

194 were used to identify and digitize the active margins of Lobe A at 1 -minute intervals (Fig. 

195 2). These areal spreading maps show that the lava lobe grew by a series of breakouts that 

196 broadly fit into two categories: narrow (0.2-0. 6 m-wide) toes that grew preferentially 

197 down-slope and broad ( 1 .4-3.5 m-wide) breakouts that formed along the sides of the 

198 lobe, nearly perpendicular to the down flow axis (Figs. 3 and 4). Initial toe thicknesses 

199 were typically ~0.2 m and, during the 35 minute period of observation, they inflated to a 

200 maximum thickness of-1 . 1 m, with most toes inflating to half their final thickness within 

201 10 minutes. Within the first minute of emplacement toes generally extended -0.2 m, with 

202 a range of 0.05 to 0.55 m. After 28 minutes, the active flow margin began to exit the field 

203 of view of the stereo cameras and toe lengths could not be accurately constrained. 

204 The mean local discharge rate measured at Lobe A over the first 19 minutes was 

205 0.0061 ± 0.0019 m 3 /s (N= 19, Fig. 5a). Note that all uncertainties in this study are 

206 reported at 1 standard deviation, a, unless otherwise noted. Standard deviation describes 
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207 how much variation there is from the mean, whereas standard error on the mean, o 

208 provides a measure of how well the mean is known based on the number of data points, 

209 N, used in the calculation. We therefore provide N as well as a to enable the calculation 

210 of er^ = a / y[N . After 35 minutes, the lobe covered a new area of 1 1.7 nr and achieved a 

a 

211 cumulative volume of 1 0.4 m . 

212 This local discharge rate was remarkably constant during the observation period, 

213 except for the interval between 19 and 24 minutes, when it decreased by nearly a factor 

214 of four. During this interval, the lobe experienced little inflation and generally lower than 

215 average areal spreading rates (Fig. 5a), which may have been due to breakouts occurring 

•5 

216 outside the field of view. At 24 minutes, the local discharge rate of 0.0061 m /s resumed, 

217 but after 28 minutes, the measured rate of volume change gradually decreased as the 

218 active flow front began to exit the field of view of the stereo-cameras. The mean areal 

219 coverage rate was 0.0055 ± 0.0038 m /s (N= 35). The areal coverage exhibits a 

220 periodicity with peaks -3 times greater than the mean (Fig. 5b) and a period of 4. 13 ± 

221 0.64 minutes (N= 8). During each cycle of areal growth (measured from trough to trough 

222 in areal coverage rates shown in Fig. 5b), the lobe covers a new area of 1.59 ± 0.28 m ( N 

223 = 7). 

224 Excess volume refers to the portion of a lava lobe that is not directly associated 

225 with its initial areal expansion. In this case, excess volume is defined as the volume 

226 change in the field of view per minute minus the change in area during that minute 

227 multiplied by the typical initial thickness of the new breakouts. The typical thickness of 

228 new breakouts from Lobe A was ~0.2 m, but may vary for other flows — generally scaling 

229 with local discharge rate (Self et ah, 1998). Excess volume changes for Lobe A appears 
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230 to be weakly anti-correlated with its areal spreading rate. However, to quantify the degree 

23 1 of anti-correlation, the cross-correlation coefficient was calculated between the areal 

232 growth and excess volume. The correlation analysis was restricted to the first 19 minutes 

233 of observation when the volume flow rate feeding the lobe was relatively constant. The 

234 areal growth rate during this 19 minute period exhibits a slight overall increase over time 

235 (slope = 0.013 nr/minute and intercept = 0.206 m 2 /minute), and therefore the trend was 

236 removed prior to analysis. After trend removal, the cross correlation coefficient, at lag = 

237 0 (the measure of direct anti-correlation) is -0.22. The negative sign of this statistic 

238 indicates anti-correlation; however, the magnitude is well within the 90% confidence 

239 interval (±0.34) indicating that the anti-correlation is very weak and not statistically 

240 significant at this level. This suggests that there may be a trade-off between areal 

24 1 spreading and inflation given an approximately constant lava supply, but that the 

242 relationship may be complicated by other factors. See Section 7.2 for further discussion. 

243 

244 6. Lava Thickness, Inflation, and the Effects of Topographic Relief on Pahoehoe 

245 Lobe Emplacement 

246 Topographic maps of the terrain before and after the emplacement of Lobe B are 

247 shown in Figure 6. Lava covered 1 1 1.83 m within the DGPS survey region, with a total 

248 lobe volume of 59.45 m . Forward Looking Infrared Radiometer (FLIR) measurements of 

249 the flow at 06:00 on 24 February 2006, revealed no anomalously hot pixels within the 

250 study region, which indicates that the flow tenninated at least 5 hours before (Harris et 

25 1 ah, 2007). Based on the time of the last field observations on 23 February, this constrains 

252 the emplacement duration to <8 hours. If the mean areal coverage and discharge rates 
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253 were similar for Lobes A and B, then the emplacement duration of Lobe B was likely 

254 between 2.6 and 5.2 hours. 

255 Using the before and after DTMs (Figs. 6a and b), topographic profiles have been 

256 extracted down-slope (Fig. 7a) and cross-slope (Fig. 7b). The lava entered the study 

257 region from the north and the lobe grew to a length of ~18 m and a maximum width of 

258 9.1 m. However, the distal part of the flow focused into a 4.3 m-wide topographic 

259 constriction before being bifurcated by a 1.0 m-high topographic obstacle into two 

260 smaller lobes that were 2.4 m and 4.0 m-wide. The profiles show that inflation of the lava 

261 tended to topographically invert the landscape, with the thickest portions of the lobe 

262 developing above topographic depressions and shallow to reverse-facing slopes, defined 

263 relative to the local flow direction. Thus, apparent barriers to the spreading of the flow 

264 are typically much smaller than the maximum height of the lobe and the thicker portions 

265 of the lobe appear confined by the thinner peripheral parts of the flow. 

266 The thickness distribution of Lobe B (Fig. 6c) above the existing surface was 

267 calculated by taking the simple difference between the final and initial topography. This 

268 yields a mean lobe thickness of 0.57 ± 0.24 m, with a maximum of 1.21 m. However, 

269 simple topographic difference overestimates inflation because lobes are initially 

270 emplaced as toes with an initial thickness. Additionally, lava had to fill small closed 

27 1 depressions before thickening above the surrounding topography. 

272 To account for the infilling of small topographic depressions, a new initial surface 

273 was created by smoothing the initial topography in a process that is analogous to tightly 

274 stretching a rubber sheet through the data. The smoothed surface was created using the 

275 GMT function SURFACE, which interpolates a grid using continuous curvature splines 


12 



276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 


in tension (Smith and Wessel, 1990). Specifically, SURFACE applies the Laplace 
transfonn (L) to the initial topography (z): 

(1-7) L x (L (z)) + 7 x L(z) = 0. [1] 

Varying the tension factor T between 1 (i.e., the harmonic surface solution) and 0 
(i.e., the minimum curvature solution) enables this data interpolation method to smooth 
the initial surface to different degrees. For instance, a grid produced with high tension (T 
— > 1) will suppress local maxima and minima by fitting an increasingly taut surface 
between the data constraints, whereas applying a looser tension factor (T— * 0) will 
produce a smooth polynomial function analogous to a flexed elastic plate that is anchored 
by the data (Smith and Wessel, 1990). 

To constrain the inflated thickness of Lobe B, a smoothed local datum was fit to 
the initial topography using a T’-value of 1 . This removes high frequency oscillations in 
the topography, while preserving the overall shape of the terrain. Subtracting this surface 
from the final topography isolates the thickness of lava above the local datum and 
provides a better proxy for inflation than the simple difference map because it reduces the 
thickness contributions due to small basin infilling. Using this method, the mean 
thickness due to inflation was 0.34 ± 0.21 m (N= 9614), with a maximum of 1.00 m (Fig. 
6d). 

To isolate local relief within the initial topography, a second interpolated surface 
was created using Eq. (1) with T= 0.25. A loose tension factor of T = 0.25 is 
recommended for low-tension interpolations to suppress extraneous oscillations (i.e., 
minima and maxima) that can otherwise be generated within the T = 0 (i.e., end-member 
minimum curvature) solution (Wessel and Smith, http://gmt.soest.hawaii.edu/ ). The new 
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299 surface preserves more variations in the initial topography than the harmonic surface 

300 solution (T= 1) and was subtracted from the original initial topography to isolate high 

301 frequency changes in elevation (i.e., relief) above and below a new local datum centered 

302 on 0 m. 

303 In the lava inundated region, the distribution of initial relief is unimodal with a 

304 mean of -0.02 ± 0.07 m (N= 9614 pixels; Fig. 8). For comparison, a 1-pixel buffer region 

305 was also defined around the outer extent of the lobe and, in this buffer region, the 

306 frequency distribution of relief had a mean of 0.00 ± 0.07 m (N= 499 pixels). However, 

307 relative to the lava inundated region, there are more pixels with a relief of 0.03 to 0.09 m 

308 and fewer pixels with a relief of -0.09 to -0.03 m (Fig. 8). This implies that there are 

309 systematically more high relief pixels bounding the pahoehoe lava lobe than in the region 

310 that was inundated. Consequently, positive relief on the order of only a few centimeters 

311 may exert an important influence on bounding the extent of low-effusion-rate pahoehoe 

312 flows. 

313 Toe-scale topography (i.e., on the order of centimeters in this case) appears to 

314 influence the areal spreading of pahoehoe lobes by blocking peripheral toes, which may 

315 then stagnate and help confine the interior portions of the flow. If the interior of the flow 

316 continues to be supplied by lava it may then inflate above the height of the surrounding 

317 barriers, thereby topographically inverting the landscape. Barriers to flow may therefore 

318 be much less than the final inflated thickness of a lobe. Based on observations of Lobe B, 

319 relief that is at least 15% of a toe’s initial thickness may be significant in terms of 

320 affecting the path of low-discharge pahoehoe lava lobes, but given sample-size- 

321 dependent restrictions on the optimum bin-size for the data (Stuges, 1926; Doane, 1976), 
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322 influences of even smaller variations in topography cannot be precluded. Obstacle 

323 heights affecting the path of higher local discharge rate flows may similarly scale with 

324 initial toe height. 

325 

326 7. Stochastic Modeling of Pahoehoe Lobe Emplacement 

327 7.1. Overview of the Model 

328 Random effects strongly influence processes of pahoehoe flow growth and pose a 

329 fundamental difficulty in developing models for their emplacement (Hon et al. 1994; 

330 Thordarson and Self, 1998; Keszthelyi et al., 1999; Crown and Baloga, 1999). Glaze and 

33 1 Baloga (2013) have recently developed a new simulation approach for modeling 

332 pahoehoe emplacement that builds on the ideas presented in Baloga and Glaze (2003). 

333 This model is based on (1) conservation of lava volume and (2) prescribed stochastic 

334 rules for lava movements within a pahoehoe lobe. This approach differs significantly 

335 from other models that treat lava emplacement as a gravity-driven viscous fluid flow on 

336 an inclined plane (e.g., Danes, 1972; Hulme, 1974; Baloga and Pieri, 1986; Baloga, 1987; 

337 Crisp and Baloga, 1990; Harris and Rowland, 2001; Rowland et al., 2004; Baloga and 

338 Glaze, 2008). The Glaze and Baloga (2013) model simulates the 3-Dimensional shape for 

339 pahoehoe lobes as they evolve in time, subject to a wide variety of ambient and internal 

340 conditions and processes. Here, the new model is used to explore the effects of inflation 

341 on pahoehoe lobes analogous to Lobes A and B. 

342 The fundamental building block of the Glaze and Baloga (2013) model is the lava 

343 “parcel”. A “parcel” of lava is defined here as a volumetric unit sampled from a 

344 probability distribution of lava volumes that typically construct a pahoehoe lobe. Thus 
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345 the smallest toe observed on the surface or along the margins of a pahoehoe lobe 

346 represents a lava parcel volume from the small end of the distribution. However, larger 

347 toes may either represent a larger parcel volume, or may be comprised of multiple 

348 parcels. 

349 A parcel is observable as a “toe” when it is permanently affixed at the surface or 

350 margin of a lobe. In the Glaze and Baloga (2013) model, lava parcels originate from a 

35 1 source region and are transported through fluid internal pathways. For the purposes of 

352 discussion here, the model assumes a constant parcel size for all parcels. Crown and 

353 Baloga (1999) measured dimensions of hundreds of toes and their mean thickness is very 

354 well constrained at -0.20 cm, (0.19 ± 0.08 m, N= 448). Based on the geometric mean toe 

355 sizes measured by Crown and Baloga (1999), the parcel volume is assumed to be 0.09 m 3 

356 (Glaze and Baloga, 2013). Assuming a square x-y grid and a typical parcel thickness of 

357 0.2 m — comparable to the observed initial thickness of toes in our study area — the 

358 corresponding cell spacing is 0.67 m, with a parcel area of 0.45 m . 

359 The Glaze and Baloga (2013) random simulation model begins with an initial 

360 source region (a single point, linear source, etc.), a constant supply of lava (i.e., one 

361 parcel per time step), and probabilistic rules governing where each parcel will be 

362 transferred and allocated. For example, in the most basic scenario of complete 

363 randomness, two random choices are made for each breakout. First, a parcel is randomly 

364 selected to be the source for the next parcel transfer. Second, a random selection is made 

365 to determine the orientation of the new parcel transfer relative to its source (i.e., North, 

366 South, East, or West). In the purely random case, every face of every parcel has an equal 

367 probability of being the source for the next parcel transfer. As various factors that 
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368 influence emplacement are examined (e.g., barriers, slope, etc.), different probabilistic 

369 rules are applied that determine where and when each parcel will be transferred. 

370 

37 1 7.2. Model Simulations of Areal Growth 

372 Figure 9a shows a typical example of a completely random simulation involving 

373 500 parcels. In this simulation, each transfer location and direction is chosen at random. 

374 Figure 10a (solid line) shows the corresponding cumulative area as a function of time 

375 step. However, Glaze and Baloga (2013) showed that topographic profdes of pahoehoe 

376 lobes are more consistent with simulations that include sequential breakouts at the 

377 margin. In this case, the momentum of new parcels breaking out from a lobe margin 

378 results in an increased probability that the next parcel will be transferred from the most 

379 recent parcel in the same direction. Figure 9b shows the influence of allowing sequential 

380 breakouts at the lobe margin on lobe morphology. Specifically, every time a new parcel is 

381 transferred to the margin (increasing the lobe area) probabilistic rules are used to 

382 determine whether 0, 1, or 2 extra parcels will be added at that location. For the 

383 simulation shown in Figure 9b, the probability of 0 additional parcels is P( 0) = 0.25, of 1 

384 additional parcel being added is P( 1) = 0.25, and of 2 additional parcels is P( 2) = 0.5. As 

385 expected, the cumulative area of the simulation that includes sequential breakouts at the 

386 margin rapidly diverges from the purely random simulation (Fig. 10a), exhibiting a 

387 steeper gradient and overall greater cumulative area. Although a general trend in the areal 

388 growth (m per time step) is evident in both simulations shown in Figure 10a, the curves 

389 are irregular with periods where the area increases rapidly and periods where the area 

390 does not change at all. Also as expected, the rate of overall areal growth tapers off over 
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391 time as the number of available lava transfer locations along the margins decreases 

392 relative to the total area of the lobe. 

393 The Glaze and Baloga (2013) simulation approach explicitly requires that lava 

394 parcels will either expand the area of the flow or contribute to inflation if they are 

395 allocated to the flow interior. This naturally results in a trade-off between areal growth 

396 and interior lobe inflation similar to the pattern observed during the emplacement of Lobe 

397 A. The natural periodicity of the areal growth rate in the simulated lobes is illustrated by 

398 the solid line in Figure 10b. Also shown in Figure 10b is the excess volume, which was 

399 defined in Section 5 as the volume in addition to what is directly associated with the 

400 change in lobe area because all breakouts will be emplaced with some initial thic kn ess. 

401 For the simulated data, the definition of excess volume results in a periodic behavior that 

402 exactly mirrors the areal growth (i.e., perfect anti-correlation). This suggests that the 

403 weak anti-correlation shown in Figure 5b may be relevant, but also indicates that other 

404 processes are present in the field that have not yet been accounted for by the model. Also 

405 in the field example, there was a period of anomalously low average lava flux between 

406 20-23 minutes, which coincided with a period of negative excess volume change. This 

407 indicates that areal spreading can also be accommodated by the depletion of lava stored 

408 within the parent lobe, rather than always being related to a simple balance between 

409 pressurization and confining strength under conditions of constant lava supply. Future 

410 studies would therefore benefit from longer time-series observations of multiple active 

411 lava lobes to further explore the statistical significance of the relationship between areal 

412 spreading and inflation. 

413 
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414 7.3. Comparisons Between Model Simulations and Observations 

4 1 5 Given the measured volume of Lobe A and typical parcel volume of 0.09 m 3 , 

416 Lobe A would include -117 parcels. The simulation of a lobe composed of 1 17 parcels 

417 begins with a 1 x 7 cell linear source region trending East-West (depicted as the gray 

418 shaded region in Figure 11a crossing through the mid-point of the array). This source 

419 region, with one parcel in each of the seven cells, is approximately equivalent to the 

420 initial length of Lobe A (i.e., -4.69 m; see Fig. 4a). The model also includes two parallel 

421 barriers located at ±3 cells along the North-South axis. The barriers are considered 

422 transitional zones in which the advancing lava enters, but becomes blocked, thereby 

423 reflecting the lava parcel back toward the interior of the flow. The distance between these 

424 barriers is 4.69 m (including the barrier cells as well), which corresponds to the 

425 topographically-bounded maximum width of Fobe A. The probabilities of sequentially 

426 adding 0, 1, or 2 extra parcels from an existing flow margin are P(0) = P(l) = 0.25, and 

427 P( 2) = 0.5. 

428 One cannot expect any specific realization of the stochastic model to exactly 

429 reproduce the morphology of an observed lava lobe, but characteristics of the simulated 

430 lobes are broadly consistent with the morphology of Lobe A. Figure lib shows examples 

43 1 of topographic profiles taken perpendicular to the flow direction through the simulated 

432 lobes averaged over 10, 20, 30, and 60 realizations. These examples demonstrate that 

433 there are stochastic variations between simulations, but also persistent morphological 

434 characteristics. For example, the simulated lobes are thickest in the middle (i.e., along the 

435 axis of the linear source region) with maximum heights that are consistent with those 

436 observed for Fobe A. The models also exhibit the effects of lava confinement along lobe 
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437 margins due to lava parcels being reflected back toward the interior of the flow. Plan 

438 views of the simulated lobes (e.g., Figs. 9a and 9b) also show irregular margins and 

439 thickness variations that are typical of inflated pahoehoe lobes. 

440 

441 8. Discussion 

442 8.1. Relationships between Pahoehoe Breakouts and Inflation 

443 Pahoehoe breakouts will be emplaced with some initial thickness and then grow 

444 through a process of inflation. This inflation process includes two dominant components: 

445 thickening of the molten core and thickening of the crust (Hon et al., 1994), with the crust 

446 consisting of both a brittle outer layer and underlying viscoelastic layer. After an initial 

447 chilled margin forms around a new toe, it will retain incoming lava and inflate as its 

448 molten core thickens to equalize the pressure with other connected portions of the fluid 

449 lava pathway. During the early stages of emplacement, the rate of core thickening will 

450 greatly exceed the rate of crustal thickening (Hon et al., 1994), but as the fluid interior 

45 1 achieves its equilibrium thickness, inflation will be increasingly dominated by the process 

452 of crustal growth (Cashman and Kauahikaua, 1997). 

453 Two conditions must be met for the upper crust to continue to accrete new 

454 material and contribute to inflation. First, there must be a continuous influx of new lava 

455 through the molten core and, second, there must be physical coupling of the molten core 

456 to the upper crust (Kauahikaua et al., 1998). If the flow of lava through the core stagnates 

457 then the lava will cool in situ and will not increase the flow’s total thickness. 

458 Additionally, if the molten core partially drains and decouples from the upper crust then 

459 new material will not be added to the base of the surface layer and the crust may even 
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460 subside and “deflate” (Kauahikaua et al., 1998). Cooling histories inferred from the 

461 thicknesses of the upper brittle and viscoelastic layers (Hon et ah, 1994) therefore 

462 provide an estimate of the time that a lava pathway flowed at full capacity, but not 

463 necessarily the total duration of emplacement if the core of the flow and the upper crust 

464 were decoupled. In some cases, the molten core of a flow can also undergo cycles of 

465 pressurization related to changes in the local discharge rate (Orr, 2011). This can lead to 

466 variations in molten core pressure that induce cycles of surface disruption and 

467 subsidence, thereby producing a range of surface textures that include slabbly pahoehoe, 

468 rubbly pahoehoe, and shatter rings. 

469 Initial cooling can strengthen the margins of a lava flow (Castruccio et al., 2013), 

470 but the process is complicated by the development of numerous cooling fractures. These 

47 1 fractures will mechanically weaken the outer brittle crust and make the underlying 

472 viscoelastic layer more important in terms of governing the overall confining strength of 

473 the flow. The frequency of breakouts will therefore depend on the balance between the 

474 combined strength of the brittle and viscoelastic layers relative to the pressure within the 

475 molten core (Hoblitt et al., 2012). If the internal pressure exceeds the confining strength 

476 of the flow’s outer layers, then the periphery of the flow will rupture and generate new 

477 breakouts that promote increased areal growth. In contrast, if the exterior of the flow does 

478 not rupture, it will pressurize and inflate through gradual deformation of the viscoelastic 

479 layer and opening of fractures in the overlying brittle crust. 

480 In this study, peripheral lobes along the margins of the PKK flow on the flanks of 

481 Klleaua Volcano, exhibited a generally constant local discharge rate (0.0061 ± 0.0019 

482 m /s). However, approximately every four minutes there was a significant increase in the 
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483 areal spreading rate. The observed periodicity in areal coverage may relate to the 

484 changing balance between the confining strength of the lava as it cools and strengthens, 

485 and the internal pressure within the core of the flow. This balance would give rise to an 

486 increased frequency of new breakouts when the core pressure exceeds the confining 

487 strength of the exterior. Enhanced areal spreading rates would then persist until the core 

488 pressure decreases below a critical threshold required to continuously rupture the active 

489 flow margin. As the rate of areal spreading deceases, continued influx of lava into the 

490 lobe would begin to re-pressurize its fluid core and perpetuate a new cycle. This process, 

491 operating on the scale of an individual lobe, is similar to what Hoblitt et al. (2012) have 

492 observed on the scale of an entire pahoehoe flow, which demonstrates the self-similarity 

493 of the inflation mechanism over a range of scales. 

494 

495 8.2. Static and Dynamic Pressure Considerations 

496 During early stages of pahoehoe emplacement, the inflation process will be 

497 dominated by molten core thickening as the liquid interior of the flow moves toward an 

498 equilibrium thickness. However, there are multiple factors that can contribute to this 

499 equilibrium. Calculations generally assume a Bingham liquid rheology, with lava 

500 descending a slope as an unconfined laminar flow under the influence of gravity (Hulme, 

501 1974; Dragoni et ah, 1986; Castruccio et al., 2013). Within this scenario, the critical 

502 depth h s required to induce down-slope flow will be 

503 h s = c 0 /pgsin(P), [2] 

504 where a 0 is the flow’s yield stress, p the flow’s density, g is gravitational acceleration, 

505 and P is the underlying slope (see Eq. [14] in Hulme, 1974). However, such calculations 
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506 of a lava flow’s equilibrium thickness will be complicated by cooling-induced changes in 

507 rheology and other factors that influence its internal pressure once the molten core is 

508 confined by a coherent crust. Several studies (e.g., Rossi and Gudmundsson, 1996; 

509 Anderson et al., 1999, 2012) have considered the “magmastatic” pressure (herein referred 

510 to as “hydrostatic” pressure) generated by fluid elevation changes along a hydrostatically 

5 1 1 connected pathway, as well as the pressure required to drive lava through pipe-like 

512 pathways and lift the crust. These models generally assume that the pathways (i.e., lava 

513 tubes) are cylindrical pipes with constant cross-sectional area, but the geometry of lava 

5 14 pathways can vary widely within a pahoehoe flow (e.g., Guest et al., 1984; Kauahikaua et 

515 al., 1998; Self et al., 1998; Calvari and Pinkerton, 1999). This raises the possibility that 

516 changes in pathway geometry can introduce additional changes in lava core pressure due 

517 to dynamic processes. 

518 These dynamic effects may be understood within the context of Bernoulli’s 

519 principle, which describes the inverse relationship between fluid flow velocity and fluid 

520 pressure within incompressible laminar flows through a pipe (Resnick and Halliday, 

521 1977; Batchelor, 1998). For lava moving within a system of combined pathways, changes 

522 in fluid core pressure could affect inflation rates by changing the magnitude of the 

523 stresses applied to the overlying crust. This may lead to increased inflation rates where 

524 the flow velocity decreases and the core pressure increases. Conversely, inflation rates 

525 would be lower where the flow velocity increases and the core pressure decreases. For 

526 otherwise equivalent lava flows through a pipe-like pathway, the most important controls 

527 on fluid velocity will be the pipe’s cross-sectional area and underlying slope. In general, 

528 lava velocities will be relatively low where the flow pathway is broad and its cross- 
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529 sectional area is large and/or where it descends a shallow slope or encounters a reverse- 

530 facing slope. Alternatively, flow velocities will be higher when lava enters a narrow 

53 1 topographic constriction and/or descends a steeper slope. These scenarios are broadly 

532 consistent with the observed relationship between the initial topography and final 

533 morphology for Lobe B. In this case, maximum inflation was observed where the flow 

534 was the broadest and above shallow to reverse-facing slopes, whereas the least inflation 

535 was observed where the flow passed through narrow constrictions and descended steep 

536 slopes (Figs. 6 and 7). These observations support the viability of a conservation of 

537 energy approach to modeling the spatial variability of inflation and its relationship to 

538 emplacement processes. They also highlight the importance of including both fluid static 

539 and dynamic components in understanding core pressure and thickening during the early 

540 stages of pahoehoe inflation. Early stage inflation processes related to core pressurization 

541 and thickening may also help to establish preferred pathways through the thickest parts of 

542 the flow, which would continue to preferentially inflate through crustal accretion if the 

543 pathways remain full. Continued crustal accretion and uplift must also be coupled with 

544 overpressurization of the molten core, which Cashman and Kauahikaua (1997) inferred 

545 from their observations of pahoehoe crustal vesicularity profiles that exhibit decreased 

546 vesicle number density with depth. 

547 

548 9. Conclusions 

549 The confining strength of a lava lobe depends strongly on the age distribution of 

550 its surfaces, with the youngest surfaces (i.e., sites of the most recent breakouts) being the 

551 weakest. Therefore, when new breakouts occur, growth tends to concentrate at these 
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552 localities. Otherwise, the next most likely location for a breakout to occur is from an 

553 older fractured surface overlying a still molten core. This process is treated with a 

554 stochastic model by introducing correlation with a probability distribution describing the 

555 preferential growth of the flow in the direction of new breakouts through sequential 

556 emplacement of lava parcels (i.e., toes) from that locality. Including correlated growth 

557 into the model also recreates observed periods of enhanced areal spreading and lobe 

558 morphologies. 

559 Time-series data for Lobe A shows that periods of enhanced areal spreading 

560 roughly alternate with periods of increased excess volume change. For an approximately 

561 constant local discharge rate, this suggests that the volume of lava not utilized during the 

562 process of areal spreading can contribute to thickening through inflation. Trade-offs 

563 between areal spreading and inflation should be detectable in thermal infrared (e.g., 

564 FLIR) imagery, provided that a single lobe is isolated in the field of view. For instance, 

565 when the flow enters a period of above average areal spreading, newly exposed lava 

566 surfaces will radiate more thermal energy and will contribute to an increased frequency 

567 distribution of high temperature pixels within the field of view. As the flow shifts toward 

568 an inflation-dominated mode, older lava surfaces will cool dominantly by conduction 

569 rather than radiation (Hon et al., 1994), and this would shift the frequency distribution of 

570 temperatures toward lower values. Temperature distributions in FLIR data may therefore 

571 be used to identify patterns of pahoehoe flow emplacement, such as periodicities in areal 

572 spreading rate and the trade-off with inflation. Statistical models of flow emplacement 

573 may be improved by coupling them with a cooling model to estimate the temperature 

574 distribution of the lava based on its emplacement age. 
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575 New toes are highly sensitive to subtle topographic irregularities and for low 

576 discharge rate flows may be deflected or stopped by obstacles measuring only a few 

577 centimeters in height. This process was observed during the emplacement of Lobe B and 

578 demonstrates that the final morphology of a lobe may inherit boundaries influenced by 

579 variations in the initial toe-scale topography that are orders of magnitude smaller than the 

580 flow’s final thic kn ess. Even though toes may be locally deflected or impeded by toe-scale 

581 topography, inflating pahoehoe lobes may overcome local barriers through a combination 

582 of stochastic breakouts in down-flow and cross-flow directions that allow them to find 

583 and follow the maximum regional slope. These behaviors are described in the model by 

584 considering the sequential emplacement of new lava parcels (i.e., toes) from existing flow 

585 margins, with boundaries that reflect parcels back towards the interior of the flow. 

586 Combining measurements of initial topography with time-series observations of 

587 active pahoehoe lobes provides new insights into the growth of pahoehoe flows. These 

588 observational constraints also help to improve statistical models of lava emplacement 

589 processes and better understand the relationship between areal expansion and inflation. 

590 This infonnation is important for understanding how subtle topographic influences can be 

591 inherited by inflated lobes to develop preferred lava pathways that influence the overall 

592 development of a flow and for interpreting the emplacement history of pahoehoe lava 

593 flows on Earth and other planetary surfaces. 

594 
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Figure 1. Lava flow map showing the study area location (19.336430°N, - 
155.055235°W) within the Prince Kuhio Kalaniana‘ole (PKK) flow (shown in yellow) on 
Kllauea Volcano, Hawai‘i. Pahoehoe lava lobes at this locality were fed lava from Pu‘u 
‘0‘o. Field measurements were acquired on 21-26 February 2006 for two lobes (A and 
B) located along the same lava pathway system. Lobe B was located directly down-slope 
of Lobe A, with a separation distance of ~10 m. The dashed lines approximately bound 
the edifices of Pu‘u ‘O‘o and Kupaianaha. Lava flow data provided by Tim Orr, United 
States Geologic Survey (USGS). 
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756 Figure 2. Orthorcctificd view of photogramme trie point cloud at time t = 35 minutes for 


757 Lobe A, with overlaid lava flow boundaries shown at 1 minute intervals and the locations 


758 of the down-slope (AA f ) and cross-slope ( BB ') profiles shown in Figure 4. 
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759 

760 Figure 3. Time series evolution of Lobe A. Dark lines show the extent of the flow at 5 

761 minute intervals with grey lines showing the flow margins at 1 minute intervals. The flow 

762 begins to move outside the field of view at 28 minutes. 
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Figure 4. (a) Down-slope (AA') and (b) cross-slope (55') profiles for Lobe A showing 


765 the changing topography of the lava lobe shown in Figures 1 and 2. 
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Figure 5. (a) Cumulative volume in the region of interest identified in Figure 2 for Lobe 
A. The mean local discharge rate (0.0061 ± 0.0019 m s" , calculated for the first 19 
minutes) is generally constant during the observation period, except for the interval 
between 19 and 24 minutes where the mean local discharge rate decreases by nearly a 
factor of four and the lobe experiences little inflation and generally lower than average 
areal spreading rates. At 28 minutes, the rate of volume change gradually decreases as the 
active flow front begins to exit the field of view of the stereo-cameras, (b) Areal coverage 
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2 

11 A rate (mean 0.33 ± 0.23 m /minute, filled black circles with solid black line) and mean 

775 excess volume change (0.23 ± 0.12 m /minute, open circles with broken line), where 

776 excess volume change is defined in the text. Excess volume changes less than zero, 
111 between 20-23 minutes, implies that areal growth was accommodated by depletion of 

778 stored lava during this time interval. The areal coverage rate shows a periodicity with 

779 peaks at intervals of 4.13 ± 0.64 minutes. All uncertainties are reported at 1 a. 
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780 

781 Figure 6. (a) Initial and (b) post-lava lobe emplacement digital terrain models (DTMs) 

782 generated from Differential Global Positioning System (DGPS) measurements for Lobe 

783 B, with cool colors showing low elevation and warm colors showing relatively higher 

784 elevation, (c) Simple difference between the final and initial topography, (d) Difference 

785 between the final topography and a Laplacian harmonic surface solution, which removes 

786 the effects of closed basin filling. Profiles CC' and DD' are shown in Figures 7a and 7b, 

787 respectively. Elevations in (a) and (b) are reported relative to the WGS84 datum, and lava 

788 thic kn esses in (c) and (d) are shown to the same scale. 
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Figure 7. (a) Down-slope and (b) cross-slope profiles showing the initial and final 


791 topography for Lobe B (see CC and DD' in Figure 6). The flow has topographically 

792 inverted the landscape, with inflation being greatest above paleo-topographic lows and 

793 anti-facing slopes, and least above forward-facing slopes. Note the pre-eruption 

794 topography and thin confining flow margins that may help to confine the thicker interior 

795 portions of the lobe. 
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796 Relief of the Pre-Emplacement Surface (m) 

797 Figure 8. Frequency distribution of relief on initial surface within the lava inundated 

798 region (dark grey) and within a 1 pixel buffer region around Lobe B (light grey). The 

799 mean relief within both regions is close to zero (i.e., 0.02 ± 0.07 m in the inundated 

800 region and 0.00 ± 0.07 m in the buffer region), but the buffer region exhibits more +0.03 

801 to +0.09 m relief surfaces and -0.09 to -0.03 m relief surface, which implies that positive 

802 relief on the order of only a few centimeters may act to confine the areal spreading of the 

803 flow. All uncertainties are reported at la. 
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Figure 9. (a) Simulation of 500 randomly emplacement lava parcels (without 
sequential breakouts from lobe margins). The maximum flow thickness is 3.2 m. (b) 
Simulation of 500 lava parcels using the same probability distribution for correlated toe 
growth used to model Lobe A (i.e., P( 0) = P( 1) = 0.25 and P( 2) = 0.5). The maximum 
flow thickness is 2.6 m. Note that including correlated breakouts in (b) increases the 
perimeter: area ratio relative to the random example shown in (a), (c and d) Plan view 
examples of two simulations of Lobe A using 117 lava parcels, boundaries (dash lines) at 
±3 cells on the North-South axis, and a probability distribution of P(0) = P( 1 ) = 0.25 and 
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813 P( 2) = 0.5 for adding 0, 1, and 2 extra parcels in the same direction when a parcel breaks 

814 out from the flow margin. The maximum height of the flows shown in (c) and (d) are 1 .2 

815 m. For visualization purposes, the REBIN function in Interactive Data Language (IDL) 

816 was used to magnify the cells by a factor of 10 and perform a default bilinear 

817 interpolation. This process smoothes the model output and provides the exterior toes with 

818 a rounded appearance. This process does not alter the raw model output shown in Figures 

819 10 or 1 lb. 
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821 Figure 10. (a) Cumulative lobe area versus time for two stochastic simulation realizations 

822 of pahoehoe lobe emplacement. The uncorrelated realization (lower solid line) was 

823 generated by considering purely random selection of both the location and direction of 

824 each parcel transfer. The correlated scenario (upper dashed line) modifies this model by 

825 introducing an increased probability of sequential parcel emplacement when new parcels 

826 break out at the periphery of the flow. In this case, the probabilities of adding 0, 1, and 2, 

827 correlated parcels are F’(O) = 0.25, .P(l) = 0.25, and P( 2) = 0.5, respectively, (b) Shows 

828 the periodicities and anti-correlation in areal spreading and excess volume change for a 

829 typical purely random simulation and an interval of 20 time steps. 

830 
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831 Distance from mid-point along North-South axis (m) 

832 Figure 9. (a) A schematic example showing how the stochastic model works. A linear 

833 source region (shown in gray) passes through the mid-point of the array and trends East- 

834 West. Two parallel barriers are shown ±3 cells from the linear source along the North- 

835 South axis. In the simulation, a cell is first chosen at random from within the confines of 


836 the flow (initially corresponding to the extent of the linear source region) and then a 

837 random cardinal direction is chosen. A lava parcel is then transferred to the specified cell 

838 (marked with a filled black circle, e.g., “1”) and to adjacent cell (identified by the 

839 associated arrow). If the parcel would be allocated beyond the confines of the flow, there 
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840 is a 0.25, 0.25, and 0.50 probability that 0, 1, or 2 extra parcels, respectively, will be 

841 allocated to sequential cells in the same direction (e.g., “2” shows 1 extra, 2 total, cells 

842 being added). As the areal extent of the flow increases, new parcel transfer locations 

843 become available (e.g., “3”). If a cell allocation would result in a parcel being extended 

844 beyond the barrier, then the parcel is reflected back towards the cell interior (e.g., “4”). 

845 Each parcel allocated to a given cell increases its thickness by 0.2 m. This process is 

846 repeated until all lava parcels are allocated, which in the case of Lobe A, includes a total 

847 of 1 17 parcels, (b) Simulated North-South profdes for Lobe A, averaged over 10, 20, 30, 

848 and 40 realizations of the model. The results show that despite the stochastic nature of the 

849 model, the average simulations converge to a stable equilibrium shape over a few tens of 

850 realizations. Also shown is the flow-perpendicular cross-section though Lobe A (Fig. 4b) 

85 1 with the underlying topographic trend removed, “i”: The maximum width simulated lobes 

852 and Lobe A are determined by the influences of topographic barriers, “ii”: Reflecting of 

853 lava parcels by the barriers produces an excess thickening in the adjacent interior cells. 

854 “iii”: High-standing topography and stagnated peripheral toes helps to confine the interior 

855 of the lobe and promote inflation, “iv”: The thickest portion of the simulated lobe is 

856 located along the medial axis and is generally consistent with the thickness of Lobe A 

857 along the flow axis. 


47 



